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ABSTRACT 

We present two large, nearby (0.013<z<0.06) samples of Early- Type Galaxies 
(ETGs): a visually classified sample of 220 ETGs, created using source- matched data 
from the Galaxy and Mass Assembly (GAMA) database with FIR/sub-mm detections 
from iJersc/ieZ- ATLAS; and a visually classified sample of 551 ETGs which are un- 
detected with if ersc/ieZ- ATLAS. Active galactic nuclei (AGN) are removed from our 
samples using optical emission line diagnostics. These samples are scrutinised to deter- 
mine characteristics of sub-mm detected versus undetected ETGs. We find similarities 
in the stellar mass distributions of the two ETG samples but testing other proper- 
ties uncovers significant differences. The sub-mm detected sample is shown to have 
lower concentration and Sersic indices than those in the undetected sample - a result 
which may be linked to the presence of dust in the former. Optical and UV-optical 
colours are also shown to be much bluer, indicating that the dust is linked with recent 
star formation. The intrinsic effective radii are on average 1.5 times larger for the 
sub-mm detected ETGs. Surface densities and groups data from the GAMA database 
are examined for the two samples, leading to the conclusion that dusty ETGs inhabit 
sparser environments than non-dusty ETGs in the nearby universe, although environ- 
ments of the brightest ETGs are shown to differ the least. Modified Planck functions 
are fit to the H-ATLAS detected PACS and SPIRE fluxes for ETGs with sub-mm 
flux densities of at least 3ct in the 350/im SPIRE band, giving a resultant mean cold 
dust temperature of Tc;=22.1K, with a range of 9-30K. The corresponding mean dust 
mass is 1.8x10''Mq, with a range of (0.08-35.0) x 10 ''Mq. The dust masses calculated 
from these fits, normalised by stellar mass, are shown to increase with decreasing stel- 
lar mass and bluer colours. Based on visual classifications of elliptical and lenticular, 
we find similar dust properties for these two early- type morphologies. We conclude 
that there is a population of elliptical galaxies which exhibit larger dust masses, lower 
Sersic index and bluer colours than the more well-known, massive, red population of 
ellipticals. 

Key words: methods: statistical - galaxies: elliptical and lenticular, cD - galaxies: 
evolution - submillinietre: galaxies 



1 INTRODUCTION 

The origin of differences between galaxy morphologies has 
aroused the interests of researchers for many decades. The 
original classifications by Hubble laid out a solid foun- 
dation for groupings of galaxies. Since then, the Hubble 
scheme has been reformed multiple times to attem pt to ex- 
plaiii characteristics associa ted with galaxies (e.g. ISandagel 
l2005l : ICappellari et al.l201ll : lKormendv fc B ender 20121). Al- 
though there is still no completely accepted revision to this 
sequence, galaxies are currently broadly classified into Early- 
Type Galaxies (ETGs) and Late- Type Galaxies (LTGs). 

ETGs form a class of typically quiescent galaxies made 
up of elliptical (E) and lenticular (SO) galaxies. The ellip- 
ticals are pure spheroidal galaxies, while SOs are spheroidal 
galaxies with a disk, both lacking spiral structure. General 
properties of this class include a tendency towards redder 
colours, little or no recent star formation, and mass and light 
distributions which decline smoothly over large radial ranges 
l|Driver et al.ll200^ ). They c an be flattened by rotation or 
due to velocity anisotropies (|lllingworthlll977h and their dy- 
namics vary, with the brighter, more massive Es shown to be 
slow rotators (EmscUcm et al. 2007). It can be exceptionally 
difficult to distinguish between Es and SOs, especially when 
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viewed face-on (|van den Berghll2009l '). hence we incorporate 
both types here. 

We intend to further our knowledge of the global char- 
acteristics of ETGs in multiple wavebands, with specific em- 
phasis on the properties and origins of ETGs which contain 
cold dust. This work is achieved by combining data from 
the Galaxy and Mass Assembly (GAMA) - a wide-area, 
multi -wavelength spectr oscopic and photometric galaxy sur- 
vey (Driver et al.ll201lh - and the Herschel Astrophysical 
Terahertz Large Area Survey (H-ATLAS) - a survey that is 
working towards providing a f ull dust census of all types of 
galaxies in the local universe (|Eales et al.ll2010l ). The com- 
bination of so many wavebands, from the FUV all the way 
through to the sub-mm 500/im band, gives a wealth of infor- 
mation that enables us to study sub-mm detected ETGs and 
their properties. This is because these wavelengths cover the 
UV/optical starlight that is absorbed by dust and re-emitted 
at MIR/FIR/sub-mm wavelengths. 

past, lack of detections in 21cm line emis- 



In the 
sion (e.g. 



. . . .Gallagher) Il972l ) and CO molecular lines (e.g. 

iBregman et al.lll992f ) in ETG spectra was assumed to in- 
dicate a lack of neutral and molecular gas and by im- 
plication, cold dust. However, dust lanes and patches 
were obse rved in optical images of some e lliptical galaxies 
(,Sadler fc Gerhardlll985l : iGoudfrooii et al.l 119941 ). Develop- 
ments in FIR and sub-mm based observations have sup- 
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ported these claims, showing that a significant fraction of 
ETGs contain cold dust within their interstellar m edium 
(ISM: lKnapp et al ll 19891 : iGoudfrooii fc de Jong||l995l ') . Some 
of t hese were subsequentl y shown to be spurious detections 
(see lBregman et al.lll99^ ). Studies at different wavelengths 
have shown this ISM to be complex: radio data led to the 
detection a nd mapping of cold n eutral gas in some ellip- 
ticals (e.g. iMorganti et al.l [20061 . and references therein); 
FIR and sub- mm data revealed signifi cant cold dust masses 
iLeeuw et al.|[20o3 : ISavov et al.ll2009l ): X-ray observations 
showed that a substantial number of ellipticals have a hot 



(« few X lO^K) plasma component (e.g. 


Kim et all 1 19921: 


iMulchaev & Jeltemal |201G|: Boroson et al. 


2011). Darticu- 


larlv in giant ellioticals PKormendv et al.l 


20091) and older 


elliDticals (Sansom et al. 2006). 



The recent influx of large-scale, multi-wavelength, 
galaxy surveys provides opportunities to delve further into 
the dust histories of ETGs. IRAS was the first infrared 
telescope to survey dust properties of ETGs, by looking 
in the MIR/FIR ban ds at 12/im, 25/im, 60nm and lOQ/im 
iBregman et al.|[l998l ). FIR emission was detected in 12-17% 
of their ETGs. Then, ISO observed deeper in the FIR , look- 
ing at ETGs at 2.5 -240^m wavelengths (|Temi et al.ll2004l : 
IXilouris e t al."2004'). The Spitzer space telescope provided 
infrared observations of some ETGs ^ mapping their distri- 



et al, in prep, Part II), we progress beyond these works with 
much larger samples of H-ATLAS detected ETGs. 

One important consideration when working with large 
samples of ETGs is determining the robustness of the sam- 
ple morphologies. How can we create a sample which is 
complete, unbiased and not contaminated by LTGs? Vi- 
sual inspection is one possibility, but has been shown to 
be subjective (jPriver et al.l |2006| ). The accuracy of eye- 
balling is also affected by the spatial and flux limita- 
tions of the galaxy photometry available from, in this case, 
the Sloan Digital Sky Survey (SDS^I). Proxies for mor- 
phology ca n be explored instead, where parame t ers such 
as co lour ([Driver et al.l l2006l: iHaines et all boOSi: iKavirail 
I2OIQI), Sersic index (heenw et al."2008': 'Bla nton et al.ll2003l : 
Driver et al. 2006: Kelvin ot al. 2012) and galaxy concentra- 
tion (|Blanton et aI.l200a : ISalim fc Richll201ol ) are associated 
with galaxy type. We investigate such proxies in the future 
paper. In this paper we make use of an eyeballed sample out 
to z<0.06 where visual classification is more reliable. 

Our investigations explore properties of sub-mm de- 
tected versus undetected ETGs. T hese include t he well- 
known morphology-density relation (|Dressleilll980l ). which 
shows that ETGs are the dominant galaxy type in high den- 
sity re gions. This h as been proven in many studies: for ex- 
ample, ICalvi et al.l (2012) recently showed that there is a 
smooth increase in the fraction of ETGs when going from 
single galaxies, to pairs, to groups. These relations may be 
explained by hierarchical merging theories, which predict 
that generations of mergers and interactions have trans- 
formed LTGs and built up to the galaxies types we see in 
the local universe. However there may be some preference 
for this build-up p rocess to occur on scales of galaxy groups 
iHovle et al.ll2012l ). 

If ETGs are formed via merging processes 
iKormendv et all |2009| ) then some may be wet or gas- 
rich mergers, resulting in residual dust in the merger 
remnants. Alternatively, could dust produced by old stars 
reproduce the dust masses currently being seen? The lack of 
correlations of cold dust mass w i th host galaxy p roperties, 
such as luminosity 



Knapdll999l : iTemi et al.lboOTi ) and age 



butions of mid- and far-IR emission (iTemi et al.l l26o^). Op- 
tical surveys include the SAURON Integral Field Survey, 
which enabled a s pectroscopic study of ionized gas in ETGs 
ijSarzi et al.ll2010l) . The ATLAS-3D volume-hmited sample 
of E TGs was used to explore their molecular gas proper- 
ties lYoung et al ] |201ll ) and showed that their molecular gas 
mass was not correlated with galaxy luminosity. A similar re- 
sult was found for neutral gas (Serra ct al. 2012). In the UV 
the GALEX mission revealed low levels of r ecent star for- 
matio n in ~ 30% of a large sample of ETGs (|Kavirai et al.l 
I2OO7I ). Additionally, extended UV disks associated with re- 
cent or potential star formation have been identified in small 
samples of ETGs (see ;Moffett et al. 2012, and references 
therein). All of these surveys showed that ETGs are not all 
quiescent, and many contain some cold dust and associated 
cool gas. 

The above mentioned FIR and sub-mm studies are now 
being superceded by the higher resolution and sensitivity 
of the Herschel Space Observatory. Herschel SPIRE and 
PACS data are being used to examine more statistically rep- 
resentative samples of ETGs. Investigations of a volume- 
limited sample of 62 ETGs, largely in the Virgo Cluster, 
with the Herschel Reference Survey (HRS), has revealed 
that some early-types contain low levels of co ol dust with 
temperatures comparable with those of LTGs (|Smith et al.l 
L2OI2). Nearby galaxies have also been examined in the 
KINGFISH project with 10 ETGs revealing the possibility 
of ongoing star formation contributing to the dust heating, 
as well as h eating from the rad iation field associated with 
older stars JSkibba et al.|[201l} ). In the local Universe, H- 
ATL AS provides the oppor tunity to work with larger sam- 
ples. [Rowlands^^e£|al] (|2012l ) used the Science Demonstration 
Phase (SDP) data to show that 42 detected ETGs (5.5% of 
luminous ETGs) co ntain as much dust as some spirals, and 
iBourne et all (|2012 f) have used stacking of GAMA sources in 
H-ATLAS sub-mm imaging to examine the properties of red 

and blue GAMA galaxies. In this paper and the next (Agius ^ lAbazaiian et al.l 1 I2OO9I) - |http: / /www.sdss.org/| 



jGeorgakakis et all 12001! ) . in ETGs older than ~2 Gyrs 
suggests that a steady accumulation of dust mass loss from 
old stars is not likely. The dust may originate externally. 
We expect most dust to be destroyed via sput tering within 
«10^ to 10* years l|Tsai fc Mathewsl 1 19951 ). if dust is 
diffusely distributed in a hot X-ray plasma. Therefore 
the detection of cold dust in ellipticals indicates that this 
process is not entirely efficient at removing all the dust 
within this timescale, or the dust is shielded within a cold 
gas phase, or is being renewed by some other mechanism. 
Thus the history of dust in ETGs is still uncertain. 

This paper utilises the GAMA I and H-ATLAS Phase 1 
catalogues, described in S'nd we assume a ACDM cosmol- 
ogy of r2Ai=0.3, r2A=0.7 and Ho=100/ikms"^Mpc"^ with 
h=G.l. We discuss the parent sample's visual classifications 
in 33 as well as additional selection criteria for the two 
samples of ETGs used in this paper: the optically selected 
GAMA/H- ATLAS matched sample and the optically se- 
lected sub-mm undetected GAMA sample. Diagnostic plots 
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and tests are presented in ^ including host galaxy proper- 
ties such as stellar mass, UV luminosities and environments. 
Dust temperature and dust mass properties are derived and 
examined in ^ Conclusions are given in 



2 GAMA AND H-ATLAS DATA 

The GAMA survey is a combined spectroscopic and multi- 
wavelength imaging campaign designed to study spatial 
structure in the nearby Universe on scales of 1 kpc to 1 
Mpc. For the spectroscopic GAMA I survey, redshifts are 
measured for galaxies in three regions of 48 deg'^ each at 9, 
12 and 14.5 hours on the celestial equator, for a combined 
observed area of 144 deg^ These regions were selected be- 
cause of the wealth of multi-wavelength data available from 
previous optical and NIR surveys. The full survey will mea- 
sure ~400,000 redshifts and their GAMA I spectroscopic 
completeness is currently 97% at SPS S Petrosian magn i- 
tude ijpet) <19.4 in all three regions (|Driver et all 1201 iT ). 
For full details of the GAMA I spectroscopic program, key 
sur vey diagnosti c s and GAMA public and team databases, 
see [Driver et al.l (|201ll '). 

The GAMA photometry campaign is built up of contri- 
butions obtained from multiple facilities and other surveys, 
including the SPS S for reprocessed optical ugriz photometry 
(York et al."2000!). GALEX for FUV and NUV photometry 
I Seib crt in prep) and the H-ATLAS survey for FIR/sub-mm 
photometry. These, in addition to other contributions, pro- 
vide moderate depth and resolution imaging data across the 
sp ectrum. The GAMA input catalogue is described in detail 
bv lBaldrv et al.l (|2010l ) and the spectroscopic tiling of these 
sources are explained in iRob otham et al. (2010). A full de- 
scriptio n of the GAMA uqrizYJHK photometry pipeline is 
given in lHill et all (120111 ). 

H-ATLAS is the widest area surve y with the Her- 
schel Space Telescope (|Eales et al.l |2010| ). To survey the 
largest possible area of sky, the maximum scan rate of 
60 arcsec sec~^ was used, with both detectors working in 
Parallel mode to obtain simultaneous observations. These 
cameras are the Photodetector Array Camera and Spec- 
trometer (PACS ), cho sen to operate at lOO/im and 160/im 
ijPoglitsch et all(2010h . and the Spectral and Photometric 
Imaging Re ceiver (SPIRE), wh ich images at 250/im, 350/im 
and 500/im l|Griffin et al.ll2010D . The angular resolutions for 
these five wavebands are approximately 9", 13", 18", 25" 
and 35" full-width half-maximum (FWHM) respectively, 
with associated 5(t point source fl ux limits of 132, 126, 32, 
36 and 45 mjy jRigbv et al.ll201ll ). 

The H-ATLAS fields were chosen to maximise the 
amount of complementary data from other surveys, which 
results in an overlap with the GAMA equatorial regions de- 
scribed above. This allows source matching between the two 
catalogues and currently gives large-scale access to ~ 10,000 
galaxies with UV, optical, NIR and FIR/sub-mm data. 

Catalogu es for H-ATLAS we re cre ated from maps a s 
described in iPascale et all (|201ll ') and llbar et all l|2010ll . 
Sources extracted from these maps had to have emission 
greater than 5a" in any of t he 3 SPIRE waveba nds, described 
in detail for the SDP in 'Rigb v et~all (120111 ). A more ex- 
tensive description of the Phase 1 dataset will be given by 
Valiante et al. (in prep). A description of the likelihood- 



ratio analysis performed to identify robust counterparts to 
the sub-mm selected so urces from the r ba nd matched cat- 
alogue can be found in ISmith et al.l (|201ll ). This uses both 
positional and photometric information of both individual 
sources and of the population in general to determine the 
reliability of an association between two sources. PACS flux 
densities are then measured by placing circular apertures 
at the SPIRE p ositions. For furthe r descriptions of these 
methods, refer to lRigbv et~al] (|201ll ). 

This paper combines GAMA UV/Optical Phase 1 
data and H-ATLAS FIR/sub-mm Phase 1 GAMA-matched 
sources to create samples of ETGs. All the internal GAMA 
data are taken from the following Data Management Units 
(DMUs) available to the GAMA team. 

• TihngCatvl6 

- CATAID, redshift quality (|Baldrv et al.ll2010l ') 

• GalexMainv02 

- FUV and NUV fluxes (|Seibertllin pr"^ 

• SersicCatv07 

- multiband Sersic derived photometry (jKelvin et al.llioi^ ) 

• StellarMassesv08 

- rest frame magnitudes, colours and stellar masses derived 
from optical SED fitting (|Tavlor et al.ll201ll ) 

• EnvironmentMeasuresvOl 

- galaxy nearest neighbour surface densities (|Brough et al.l 
submitted^ 

• GamaGalaxyGroupsCatv04 

- galaxy groups and multiplicities (|Robotham et aLll201ll ) 

• kCorrectionsv02 

- GAMA-I multi-band k-corrections 



3 EARLY-TYPE GALAXY SAMPLE 
SELECTION 

To create our nearby sample of ETGs, we took a volume- 
limited, visua l ly classi fied sample of GAMA-I galaxies from 
I Kelvin et al.l (|in pred ). This is restricted to the redshift 
range 0.013<z<0.06, a Galactic extinction-corrected limit 
of Tpet <19.8 and an additional absolute magnitude cut at 
Mr <-17.4 to ensure a volume-limited sample. After an ad- 
ditional 19 galaxies were removed due to poor imaging data 
or bad SDSS photometry, this resulted in a parent dataset 
of 4,110 galaxies, with a range of galaxy types ([Kelvin et al] 
lin pred). 

3.1 Visual Classification 

The visual classification for this parent sample wa s carried 
out and is fully described in Kelvin et al.l (|in prepi) . Briefly, 
they generate three colour postage stamps using UKIDSS 
H and SDSS i and g bands for red, green and blue colours 
respectively. These were then classified on a step-by-step 
basis by three independent people across three tiers. The 
tiers were defined as follows: 

Is the galaxy spheroid- or disk-dominated? 
Is it a single- or multi-component system? 
Is one component a bar? 

These criteria allowed an analysis of the broad-scale mor- 
phological properties of the galaxies, independent of colour. 
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0.1 1.0 10.0 -1.5 -1.0 -0.5 0.0 0.5 

Sersic Index (r band) log„([NII] ?L6583/Ha) 



Figure 1. Optical colour and Sersic index properties of the visu- 
ally classified parent sample, coloured by galaxy type. 



Classifications were assigned in cases where at least two of 
the observers agreed, thereby allowing for objective visual 
classifications. Spheroid-dominated systems are all classi- 
fied as ETGs, and then further subdivided into E (single- 
component) galaxies, SOa (multiple-component) galaxies 
and SBOa (multiple-component with a bar) galaxies. Note 
that the SOa grouping includes all sub-catagories of SO, SO/a 
and Sa galaxies. Spiral galaxy classifications, which are not 
the subject of this work, are assigned to systems which are 
classed as disk-dominated. For example, Sd galaxies are clas- 
sified as disk-dominated, single-component systems. 

The distribution of this parent sample in optical colour 
and Sersic index space is shown in Fig. [T] where Sersic 
indices are extract ed from GAMA Sersic photometry fits 
l|Kelvin et al.ll2012l ) and colours are the restframe colours 
taken from Stellar Population Synthesi s modelling of opt i- 
cal spectral energy distributions (SEDs: iTavlor et al.ll2011f ). 
The classifications resulting from the previously described 
eyeballing are illustrated by colours identified in the colour 
bar. This plot gives us an indication of the properties of 
each galaxy type. It shows that as we go to later types, 
both colour and Sersic index properties decrease on av- 
erage. This can be interpreted as earlier types tending 
towards redder colours and less exponential Sersic pro- 
files (n>l), a r esult co nsistent w ith previous stud ies such 
(I2OO2I ): 



Peng et al 



as IFeng 
(l2Q0d);|l 



Blanton et a l. (2003); Dr iver et al.l 



Haines et all l|2008l ): iKavirai et al.l (|201lh . 



3.2 Removal of Active Galactic Nuclei 

As the key focus of this work is to compare the properties 
of two samples of ETGs, it is necessary to make sure the 
properties obtained accurately represent the host galaxies. 
In addition, we are not studying the dust heating contribu- 
tion, or sub-mm emission, from AGN. For these reasons we 
now discuss our AGN s election method and removal. 

[Foster et al] bOl'j) use the Gas and Absorption Line 
Fitting (GANDALF: ISarzi et al1l2 006l') algorithm to measure 
stellar emission lines from fiux calibrated GAMA spectra. 
They define galaxies with significant emission lines as those 



Figure 2. This BPT diagram shows the flux ratios of [OIII] 
A5007/H/3 emission lines against [Nil] A6584/Ha emission for 
all GAMA emission line galaxies (grey points), with the 
iKaufl^mann et al. AGN dividi ng line shown as the solid 

black line and the lKewlev et al] 1I2OOII ) dividing line as the dashed 
black line. Visually classified galaxies with CATAID-matched 
emission lines are shown and coloured according to their mor- 
phological classification. 



with 3cr levels in Hq, H/3 and [NII]A6 584 lines. With their 
results, we produce a BPT diagram (|Baldwin et aLlfigSlI ) 
for GAMA galaxies with these lines (see Fig. |2]). Emission- 
line galaxies in our parent sample are over-plotted to indi- 
cate their distribution and are again colo ured by eyeballed 
morph ology. Based on the prescription bv lKauffmann et alj 
(120031 ) (Eq. 1), we show the dividing line between star- 
forming galaxies and AGN-dominated galaxies, in the form: 



log 



[OIII]\ 
HP ) 



0.6 



lof 



+ 1.3 



0.05 



(1) 



This line identifies the upper limit for emission line mea- 
surements which are representative of pure star formatio n in 
the BPT plot. Fig.[2]also shows the lKewlev et al.l (|200ll ) up- 
per limit line (Eq.[2]) for extreme starburst limits. 



loe 



([OIII] 



HP 



0.61 



log 



-I- 1.19 



0.47 



(2) 



Note that this latter relationship is not our AGN identifier 
and is only shown for comparison as a soft limit. All galaxies 
above and to the right of the Kauffmann dividing line are 
considered to have optical AGN emission. 

When matching our parent sample to the GANDALF 
emission lines catalogue, we get 1651 (40%) matches for 
emission lines. From this selection, 343 (21%) galaxies are 
dominated by AGN and 1308 (79%) are star-forming. Lack 
of significant emission lines is the chief cause of unmatched 
galaxies. The colour selection in Fig. [2] indicates how these 
are divided by morphology and we show the numerical re- 
sults in Table [T] The first row in Table [T] shows that early- 
types (E, SOa and SBOa) have larger proportions of their 
emission-line galaxies identified as AGN (~59%) than later- 
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E 




SOa 




SBOa 




Sbc 


SBbcd 




Sd 




SF 


AGN 


SF 


AGN 


SF 


AGN 


SF 


AGN 


SF 


AGN 




AGN 


Emission-line galaxies 


71 


83 


20 


18 


1 


3 


74 


10 


6 


5 


1136 


224 


H- ATLAS detected 


2 


2 


7 


6 





1 


39 


5 


6 


4 


27 


3 


H- ATLAS undetected 


69 


81 


13 


12 


1 


2 


35 


5 





1 


1109 


221 



Table 1. Galaxies in our parent sample with emission lines are divided into star- forming (SF) and active galactic nuclei emitting (AGN). 
Numbers shown here are divided into H-ATLAS 250^m flux detected (102 galaxies) and those which are undetected by H-ATLAS (1549 
galaxies). They are additionally separated into their morphological classifications. 



types (^^26%). However, both Table[T]and Fig. [2]show that 
late-types preferentially occupy the BPT diagram and for 
our parent sample, emission from these late-types is most 
likely powered by young stars rather than AGN. 

For the next stage of our sampling, we remove all galax- 
ies with optically identified AGN emission from our parent 
sample, leaving us with 3767 galaxies to pick ETGs from. 
The following section d escribes how we use t hose galaxies 
classified as E or SOa by I Kelvin et al.l (|in prepl) to form two 
samples which are the subject of this work. 



3.3 H-ATLAS Detected and Undetected Samples 

We chose to create two samples out of the remaining vi- 
sually classified ETGs: the optically selected, sub-mm de- 
tected, ETG Sample (SubS), and the optically selected, sub- 
mm undetected, ETG Sample (OptS). SubS is the eyeballed 
sample of E (elliptical) and SOa (lenticular/early Sa spiral) 
galaxies which have positional matched detections to within 
1" in both GAMA and H-ATLAS. OptS is the eyeballed 
sample of GAMA ETGs which do not have H-ATLAS de- 
tections (5(T in any SPIRE waveband). As described, our 
samples contain Sa galaxies, but we make the assumption 
that they are not numerous and will not significantly skew 
our results. This is likely due to additional imposed selec- 
tion criteria (described below), together with published Sa 
popu lation percentages (~18% of ETGs; iNair fc AbrahamI 
I2OIOI . their Table 3). 

We also imposed additional selection criteria to form 
these samples. This includes setting an ellipticity cut of 
(1 — (6/ffl)) < 0.7 to remove any edge-on disk galaxies. Addi- 
tionally we ensured that all our sample ETGs had an efTec- 
tive radius of at least the seeing FWHM (1.2") to minimise 
any uncertainty in parameters related to seeing. These cri- 
teria removed 32 and 62 galaxies from the SubS and OptS, 
respectively. As described in i\3.2l we removed galaxies char- 
acterised by optical AGN emission, showing these numbers 
in rows 2 and 3 of Table [T] for H-ATLAS detected and un- 
detected galaxies respectively. For the H-ATLAS detections, 
we also removed 15 galaxies which had low reliabilities and 
are therefore considered likely to be false counterparts. Fi- 
nally, we ran an independent visual check on these galaxies 
to remove any with obvious spiral arms. This was a neces- 
sary additional criteri on as spiral arms were not accounted 
for in the eyeballing bv lKelvin et al.l l|in prepl) . and 22 galax- 
ies were removed in this process. Our resultant samples in- 
cluded 229 ETGs (33% E and 67% SOa) in the SubS and 
551 ETGs (54% E and 46% SOa) in the OptS. Therefore the 





Figure 3. Example images of sub-mm detected galaxies with E 
classifications (top row) and SOa classifications (bottom row) . The 
images are 20" SDSS g band cutouts. 



sub-mm detected ETGs are more dominated by lenticulars. 
Examples of SubS ETGs are shown in Fig. |3l 



3.4 Selection Effects and Completeness 

Both SubS and OptS are affected by selection effects and 
completeness issues. In this section we will examine these 
and try to assess how they will affect our comparisons. 

First of all, both samples are affected by a lack of 
GAMA catalogue ID (CATAID) completeness. The 229 and 
551 ETGs in the SubS and OptS respectively all have 
GAMA CATAIDs and redshifts with normalised redshift 
qualities nQ>3 (good for science). However, due to the fail- 
ure in extraction of derived parameters for a small subset of 
galaxies in the GAMA database, our samples do not 100% 
match CATAIDs in all GAMA internal catalogues used in 
this work. In particular, NUV magnitudes are not complete 
for both samples. We explore the level of this incompleteness 
in UJ\ 

Fig. |4] shows the GAMA absolute r band magnitude of 
the two samples, as a function of redshift. The plot marks 
the r band spectroscopic complete ness limit of tpet = 19.4 in 
all three fields as the dashed line ([Driver et al ] l201l[ ). Note 
that there are galaxies plotted below this limit which will be 
as faint as rpet = 19.8: the magnitude limit for the ongoing 
GAMA-II data collection. This plot also highlights the ab- 
solute magnitude cutoff used here to form a volume-limited 
sample (down to M,. = — 17.4mag., i.e. slightly brighter than 
the Small Magellanic Cloud) . Less luminous galaxies are re- 
moved from our samples in this paper. A volume-limited 
sample going out to higher redshift would require a brighter 
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Figure 5. The orange histograms show (a) 250/im, (b) 350/im and (c) 500/im H-ATLAS flux values for our sample of optically-selected, 
early-type gal axies with 5cr det ections in at least one band. The short dashed line gives the 3cr flux level, and the long dashed line the 
5(T flux level l lEales et al.ll2010l '). The distribution of the SubS is above these levels in the 250/im regime, due to the selection criteria, 
but '~35% of 350/im fluxes are below the 5cr level and most of the distribution lies below the 5cr flux level for the 500;im band. 
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Figure 4. The distribution of all GAMA detections (grey points) 
plotted as a function of redshift up to z=0.2. The dashed line 
shows the GAMA r band spectroscopic completeness limit of 
mr=19.4. The vertical dashed lines illustrate the redshift cut- 
offs for our early-type samples and the horizontal dashed line 
indicates the absolute magnitude cutoff. The sub-mm detected 
sample is shown as blue diamonds and the undetected sample as 
black diamonds. 



magnitude cut, an issue which will be dealt with in Paper 
II of this series. 

This plot shows an apparent lack of completeness at 
higher redshift for our samples. The parent sample is com- 
plete even in this regime and therefore this effect is a result 
of the methods used to create the two samples. We investi- 
gate this phenomenon on a galaxy-by-galaxy basis and find 
the removal of AGN and galaxies with low effective radius 
are the main cause of this incompleteness. However, this ef- 
fect is true for both samples, indicating that comparative 
measures between the two are unlikely to be biased. 

We now discuss completeness issues specific to SubS. 
These sources were selected based on H-ATLAS SPIRE de- 
tections greater than 5a in any waveband. As shown in Fig. 
[5ja), all these 5cr detections are in the 250/im waveband. A 



study in lRigbv et all (120111 ) shows H-ATLAS SDP as having 
greater than 80% catalogue number density completeness. 
This missing ~20% is due to a number of undetected faint 
sources indicated by random noise fluctuations in the sim- 
ulated maps, or because of source blending. This is likely 
to be larger f or the Phase 1 data due to its deeper region. 
iDunne et al.l (|201ll ) also shows a slight breakdown in ID 
completeness, as they cannot guarantee that all detections 
which are given SDSS ID counterparts have been correctly 
identified. This is due to positional uncertainties, close sec- 
ondaries and the random probability of finding a background 

source within 10". 

Based on recommendations by ISmith et al.l (|201ll ) , all 
of our sub-mm galaxies have a reliability R>0.8, where reli- 
ability is the likelihood that each object is the correct coun- 
terpart out of all the counterparts within a search radius. 
This ensures that the contamination rate is kept to a mini- 
mum and that the SDSS r-band source is linked to the FIR 
emission. Based on this value, we are able to estimate the 
likely numb er of false IDs in t he SubS. This is done using 
Eq. 12 from lSmith et all l|201ll) : 

Nifalse) = ^ (1 - J?) (3) 

B>0.8 

This results in ~3 galaxies (1%) in our sample having a 
possibly false ID. We judge this fraction low enough to not 
be a major problem for the rest of this work. 

We test the galaxies in SubS to check for the likelihood 
of contamination from strong lensing sources. We use two 
methods to ch e ck for lenses. The first method is a study by 
iNegrello et"al1 l|2010l ). where galaxies with SPIRE 500^m 
emission > 100 mjy are likely to be strong lenses. Seven 
of the galaxies in the SubS fulfill this crite r ion. T he sec- 
ond method is that o f iGon zalez-Nuevo et al] (|2012l ). where 
galaxies fulfilling all the following criteria have a 50% likeli- 
hood of being strong lenses. 

SPIRE 250 ^m fiux emission > 35mJy 

SPIRE 350 yitm fiux emission > 85mJy 
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SPIRE 350/250 ^^m ratio > 0.6 
SPIRE 500/350 ^im ratio > 0.4 

Based on this second set of criteria, we find two ETGs in the 
SuhS are 50% likely to be strong lenses. However, these two 
methods do not pick out the same galaxies as being possible 
strong lenses. Altogether, these tests indicate low (^^4%) 
levels of contamination from lensing, but we remove these 
nine galaxies from the SubS as we want to avoid including 
any possible strong lenses or background contaminants in 
our sample. This reduces our SuhS to 220 ETGs. 

The second and third plots in Fig. [S] indicate that in 
the 350/im band, 65% (142 detections) of the galaxies are 
above 5cr and only 14% (31 detections) are above this limit 
in the 500/im waveband. This causes some problems, as will 
be discussed later in ij5] with fitting SEDs to the data to 
obtain dust-related parameters. 

In subsequent sections everything is examined on a com- 
parative basis between the purely optical and optical-f-sub- 
mm samples whenever possible. This comparative approach 
allows us to investigate the relative behaviours of sub-mm 
detected versus undetected ETGs in the nearby universe. 



4 SUB-MM DETECTED VS UNDETECTED 
DIAGNOSTICS 

iKormendv et al.l (|2009^ split elliptical galaxies into two main 
classes: Giant Ellipticals and Normal/Dwarf Ellipticals. The 
primary class can be defined as having a brightness M,- < 
—21.5 and n>4, whilst the Normal/Dwarf class has the re- 
verse criteria. To gain understanding of the characteristics 
of galaxies seen in our samples, we plot distributions of mul- 
tiple parameters in Fig. |6] for both our Sub-mm and Optical 
samples. At a glance, most of these distributions ar e seen to 
differ for the SubS and OptS. Based on the Kormendv et al.l 
brightness definition, both samples contain a broad 
range of r band absolute magnitudes, indicating a mix of 
both giant and normal ETGs. 

To compare the structural and intrinsic properties of 
ETGs in our samples, we carry out a series of Kolmogorov- 
Smirnov (KS) tests for our samples and subsets thereof, the 
results of which are shown in Table [S] We choose to test 
the samples as a whole because they are volume-limited, 
and therefore the testing is unlikely to be biased by, for ex- 
ample, redshift or mass effects. We also test two subsets: 
the first containing those galaxies classed as Giant ETGs 
(Mr <-21.5) and the second classed as Normal/Dwarf ETGs 
(Mr >-21.5). The aim of this subset testing is to investigate 
whether these two so-called separate populations have their 
own unique set of properties for sub-mm detected and unde- 
tected ETGs. We find that 36% of Giant ETGs are sub-mm 
detected, whereas only a quarter of the Normal ETGs are 
detected. 

Table [5] therefore shows parameters which appear to be 
most representative of the characteristics of the two samples 
of galaxies. We use a combination of Table [2] and Fig. [6] to 
investigate how these distributions vary for our two samples. 
We choose the KS-test because of its usefulness as a non- 
parametric test for checking a null hypothesis for unbinned 
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Absolute Magnitude (r band) 

Figure 7. Scatter plot showing the colour-magnitude (CM) dis- 
tribution of galaxies in the SubS (blue circles) and OptS (black- 
grey circles). The division between Red Sequence and Blue Cloud 
defined in the text is shown as the black solid line. 



distri butions tha t are functions of a single independent vari- 
able ()Presslll99j ). We choose KS-probabilities of 1% as the 
significance level for our distribution testing. 

By examining the spectroscopic redshift values in this 
table, we can deduce that the typical redshift for the galaxies 
in both samples do not differ for our whole samples or Giant 
subset, although there is a difference in distributions for the 
Normal subset. 

Concentrati on index, defined a s the ratio of Petrosian 
r90 to r50 radu (|Mateus et al.lbOOd ). where r90 and r50 are 
the SDSS circular aperture radii within wh ich 90% and 50% 
of the flux are contained (|Bell et al.ll2003l ). is the first pa- 
rameter plotted in Fig. [6l It is linked to the concentration 
of light in the centre of a galaxy, and it ha s higher values i n 
ETGs than in LTGs (|Strateva et al.ll200ll : [Conse hcc"2006l'). 
Both Table [2] and Fig. [HJa) show that concentration index 
for galaxies in the two full samples differs significantly, with 
the probability of getting the null hypothesis of no difference 
approaching zero. The distribution for the SubS is offset to- 
wards low concentrations for the full samples and for both 
low and high luminosity subsets. These differences indicate 
that H- ATLAS detected ETGs are likely to have less concen- 
trated light distributions, and therefore may have more per- 
turbed internal structure than those which are undetected, 
indicating some past merging or formation activity. 

We next l ook at Sersic profile properties for the samples. 
iKelvin et~al] l|2012t ) describe how their SIGM A wrapper 
around the GALFIT3 profile fitting program (IPeng et al.l 
i2010 ) is used to fit the entire GAMA database with 2D Sersic 
profiles in multiple wavebands. The key parameters that we 
are interested in from this profile fitting are the Sersic in- 
dices and effective radii of galaxies in our samples. Initially 
we examine Sersic index, which gives us information on the 
distribution of light within a galaxy. This parameter is more 
useful than concentration index as it is less affected by see- 
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Figure 6. Histograms representing distributions of ETGs in the SubS and OptS for the optical parameters discussed in §|4] From left 
to right, top to bottom, we show the distributions of concentration index, Sersic index, intrinsic effective radius, (u-r) colour, absolute 
r band magnitude, and stellar mass. Plots include average errors for these parameters, as well as error bars for the normalised binning. 
KS- probabilities from Table [2] are also included for each set of distributions. 



ing, having been convolved with the point-spread function 
(PSF), and therefore does not vary so much with redshift. 
The recovered Sersic index distributions show similar trends 
to those observed for concentration index. The distribution 
of the samples in equal log spacing shown in Fig. [6{b) indi- 
cates that ETGs in the SubS have lower Sersic indices, and 
KS-testing the samples shows them to be significantly dif- 
ferent. One possible interpretation of this result is that the 
difference found in the two sa mples is due to the effect of 
dust, which has been shown (eg. lPastrav et al.l2013l ) to lower 
values of measured Sersic index. Alternatively, these results 
for both Sersic and concentration index may be caused by 
morphological disturbances, or problems with the morpho- 
logical selection process, which is limited by the resolution 
of the SDSS. 

It is interesting to note that mean Sersic indices for the 
two samples are higher for the Giant ETG s subset than the 
Normal ETGs subset, as is predicted by the lKormendv etld] 
(|2009l ) separation criteria. In addition, the observed differ- 
ences in distributions of Sersic index affects mainly the Nor- 
mal ETG subset (see Table [2]). Thus, if we believe that this 
difference is caused by dust, Sersic index would appear to be 



most affected by its presence in low luminosity ETGs. This 
could indicate that higher luminosity ETGs contain lower 
normalised dust masses, leading to higher Sersic indices and 
profiles more similar to OptS ETGs. We explore this effect 
with respect to specific dust mass in ^ 

Apparent effective radii are also calculated for the 
galaxies through this Sersic profile fitting, and we convert 
these to intrinsic radii using simple geometry and angu- 
lar diameter distances. The distributions are plotted in log 
space in Fig.[SJc). It appears that OptS ETGs typically have 
smaller effective radii than those in the SubS, although this 
result is not robust when the sample is divided into lumi- 
nosity subets. 

We examine (u-r) colours of both samples, where galax- 
ies in the optical Red Sequence (RS) are separated from the 
Blue Cloud (BC) using the following relation: 

u-r = -0.09Af,. - 0.0347 (4) 

This was calculated by fitting approximate straight lines to 
the GAMA RS and BC, and finding the line equidistant in 
magnitude from these fits. This is a similar method to other 
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Table 2. Two-tailed Kolmogorov-Smirnov (KS) test results for the distribution of host galaxy parameters for our ETGs, tested as 
follows: (1) testing the full samples, (2) testing those ETGs in the samples with Mr <-21.5 and (3) those ETGs with Mr >-21.5. The 
first column gives the parameter which has been tested, with the two distributions being the sub-mm detected and undetected samples 
respectively. The second column gives the KS-statistic D, which is defined as the maximum value of the absolute difference between 
two cumulative distribution functions. The third column gives KS-probability, which is the probability for the null hypothesis that these 
data sets are drawn from the same distribution. The fourth and fifth columns show the numbers of SubS and OptS galaxies tested. The 
final two columns show the mean values for the SubS and OptS parameter distributions. The test is carried out using the IDL routine 
KSTWO. 



publ ished works for different optical colours (e.g. iBell et al] 
I2OO3I V ETG colours are explored in Fig. [71 where the (u- 
r) CM scatter plot is shown for both samples. The diagram 
highlights the point that ETGs, whether detected in the sub- 
mm or not, have a multitude of both blue and red colours. 
Interestingly, the SubS does not show as dense a cluster- 
ing in the Red Sequence as the OptS and is more evenly 
distributed towards the blue end. Both Figs. [6l^d) and [7| in- 
dicate that the SubS contains a larger proportion of blue 
ETGs. Table [2] also shows the distribution of optical colours 
is quite different, even when testing separate subsets. 

The full sample testing of absolute magnitudes indicates 
different distributions. However, results from subset testing 
indicate that the Giant ETGs have similar distributions, 
whilst Normal ETGs are different with ~0.6% probability. 
This is unsurprising, given the OptS ETGs extend to fainter 
magnitudes than SubS ETGs, shown in Fig. [71 

Given the link between luminosity and stellar mass of 
galaxies, we would expect equivalent results for stellar mass 



distributions. However, all the samples have similar mass 
distributions, although once again the result is less robust for 
the Giant subset. On average, the SubS reaches to brighter 
absolute magnitudes than the OptS, whereas mass differ- 
ences between the samples are not significant. This may 
mean a decrease in mass-to-light ratios for SubS ETGs. 

Altogether, after testing these intrinsic properties, we 
find sub-mm detected E TGs to be bluer tha n undetected 
ETGs (see also results bylPariu sh et al1l201ll '). and we also 
show they are likely to have similar masses, lower Sersic in- 
dices, are less centrally concentrated and less compact than 
undetected ETGs in our samples. 

4.1 UV Parameters 

Availability of GALEX data means that we can look at the 
ultraviolet (UV) bands for both our OptS and SubS, typ- 
ically used as an indicator of recent star formation. The 
GALEX bands are Far UV (FUV), with an effective wave- 
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Red Sequence 


50 0.27 


373 0.78 
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88 0.48 


77 0.16 


Blue Cloud 


46 0.25 


31 0.06 



0.20 



Table 3. Numbers and Fractions of ETGs residing in the UV- 
O ptical Red Sequence, Green Valley and Blue Cloud, as defined 
bv lBourne et al l ||2012| '). 



length of 1528A, and Near UV (NUV), with an effective 
wavelength of 2271 A (|Martin et al.l 12005'). These GAMA- 
GALEX catalogues seek to reconstruct the original UV flux 
for each given optical source, and include UV fluxes and AB 
apparent magnitudes for both wavebands. 

We corrected NUV CATAID-matched apparent mag- 
nitudes for Galac tic extinction (for more details, see 
IWvder et all bOOTl l and converted them to absolute mag- 
nitudes using GAMA database k-correc tions derived from 
kcorre ct v. 4- 2 software, described in iBlanton &: Roweii 
|20o3). This IDL code calculates K-corrections for the galax- 
ies based on the best flt sum of templates to an SED and 
accounts for both the angular diameter distance and cosmo- 
logical surface-brightness dimming for the distance modulus. 

We obtain 481 and 184 galaxies with GALEX NUV de- 
tections in the OptS (87%) and SubS (83%) respectively. 
The SubS is brighter on average in the NUV by ~2 mags. 
We examine the NUV-r colours in Fig. [S] where we see the 
SubS contains bluer UV-Optical colours than the OptS, sim- 
ilar to the H-ATLAS resuh for ETGs in the SDP field of 
[Rowlands et all (|2012f ). As in ^ we perform KS-tests on 
the full SubS and OptS samples and their luminosity sub- 
sets. The results of these tests are presented in Table [5] and 
show the distributions of UV colours are very different for 
ETGs in all these cases, with SubS ETGs showing much 
bluer colours overall. 

For comparison, we also show where H-ATLAS detected 
LTGs (visually classified in the same process as ETGs in our 
samples) lie on this histogram. The colours of these LTGs 
are bluer on average than both our ETG samples, but the 
SubS straddles the gap between blue and red modes filled by 
LTGs and OptS respectively. This indicates that H-ATLAS 
detected ETGs are forming a colour population of their own. 

We next examine the UV-optical CM diag ram in Fig.|9l 
and follow the method of iBourne et all (|2012l ). using NUV- 
r colour boundary cuts shown in Eqs. O [H] and [7] below to 
define regions as the Red Sequence (NUV-RS), Green Valley 
(NUV-GV) and Blue Cloud (NUV-BC) respectively. 

1.23 - 0.17M,. < (NUV ~ r)rest < 7.0 (5) 
-0.27 - 0.17M,. < {NUV ~ r)rc^st < 1.23 - 0.17Mr (6) 
1.0 < {NUV - r)rest < -0.27 - 0.17Af,(7) 

The results of applying these boundary conditions to de- 
fine different regions are shown in Table [3] This highlights a 
key difference between the SubS and OptS. There is a large 
fractional difference between the populations in all these 
regimes, with OptS galaxies most prominently located in the 
NUV-RS, and SubS ETGs mostly dominated by the NUV- 
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Figure 8. Histograms of UV-Optical colour for OptS (black) and 
SubS (blue), as well as H-ATLAS detected LTGs (grey). The SubS 
appears to be skewed towards the bluer end, whilst the peak of 
the Opts is decidedly within the red end of the plot. The blue, 
black and grey dot-dashed lines show the respective mean colours 
for the SubS, OptS and LTGs (4.05, 5.33 and 3.31 mags). 
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Figure 9. The distribution of SubS (blue filled circles) and OptS 
(black filled circles) on the UV-optical colour -magnitude diagram . 
Dashed lines indicate the cutoffs taken from lBourne et alj 1I2OI2I I 
to separate galaxies into Red Sequence, Green Valley and Blue 
Cloud respectively. 



GV. It is difficult to determine whether the 48% sub-mm 
detected galaxies in the NUV-GV are ETGs in transition 
from the Blue Cloud to the Red Sequence, or whether they 
are NUV-BC ETGs that are dust-reddened to appear in the 
NUV-GV, although absolute NUV flux levels do suggest on- 
going or recent star-formation in SubS ETGs. 

The sub-mm detected ETG occupation of the NUV-BC 
is a factor of three more than that of undetected ETGs, pos- 
sibly indicating more ongoing or recent star formation in the 
former. If the UV emission in ETGs is dominated by emis- 
sion from recent star formation and dust is predominantly 
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heated by photons from the young stellar population, then 
we would expect SubS ETGs to be dust-reddened and driven 
towards the NUV-RS. The fact that the SubS ETGs occupy 
mostly the NUV-GV implies that one of these assumptions 
may not hold. Additionally, we might predict that the dust 
does not lie within the same geometry as the star formation, 
potentially implying star formation is not the main source 
of heating the dust. 

4.2 Environment Parameters 

The iDresslei] (|l98Cll ) morphology-density relation is a key 
driver for investigating the varying environments of ETGs. 
Although it has been shown that bulge-dominated galaxies 
are more commonly found in the densest regions of the Uni- 
verse, it is unclear whether ETGs with differing dust prop- 
erties are typically found within different environments. 

Disk-dominated galaxies are thought to be transformed 
into bulge-dominated galaxi es via processes which o ccur 
in dense environments (see iBoselli fc Gavazzil (|2006h for 
a review of such processes). This work is only concerned 
with bulge-dominated galaxies, or ETGs. However, we are 
comparing sub-mm non-detected passive ETGs, which may 
likely be at the end stages of their lives (OptS) with dusty 
ETGs which are likely to have active star-formation (SubS). 
It would be interesting to see whether these two classes of 
ETGs reside in different environments. 

The surface density information used to investigate 
environments for our ETG samples was taken from the 
GAMA EnvironmentMeasures database. These surface den- 
sities (Sgoi) in galaxies Mpc~^ were calculated based on the 
Nth Nearest Neighbour method, using: 

where N is the nearest neighbour number, set at a 
value of 5, and Djv is the projected distance in co- 
moving Mpc to the N"" nearest neighbour within a veloc- 
ity cylinder ±1000kms~^ from a volume-limited, density- 
defining population. This population is defined by an r 
band absolute magni tude Mr <-20 and a redshift selec- 
tion of 0.0 02<z<0.18 (jWiiesinghe et al.ll2012l : lBrough etld] 
Isubmittedl) . We choose to include upper limits calculated 
based on the distance to the nearest angular edge. Only one 
galaxy within our two samples is removed because the den- 
sity cannot be estimated, reducing the SubS to 219 ETGs. 
We find surface density to be an extremely useful parameter 
for us to work with because it gives us information about 
the area around the galaxy and whether there are likely to 
be any real interactions occurring between it and its neigh- 
bours. 

Fig. 1101 shows the surface densities for both the reduced 
SubS and OptS. The OptS is more distributed over the range 
of surface densities, indicating more of these galaxies are 
in denser regions. As in ^ we KS-test the distributions 
of densities for the samples and within luminosity subsets. 
Results given in Table [2] show that the ETGs have different 
density distributions in the main sample, as well as both 
luminosity s ubsets, w it h pro babilities of below 1%. 

Robotham et al.l (|201ll ') produced a GAMA Galaxy 
Group Catalogue (G'^C) that identifies which groups the 




Surface Density (Galaxies Mpc") 

Figure 10. Normalised histograms of surface density in log bins, 
where the blue histogram is the SubS and black histogram is the 
Opts. The error bars represent average errors on both surface 
density and the normalised binning. KS-probability for compari- 
son of surface density distributions from Table [2] is also included. 



GAMA galaxies belong to, as well as assigning each group 
a set of properties. It is assumed that if the galaxies are not 
assigned a Group ID, then they are in the field. We have 
used G^C to assign our galaxies either field or group sta- 
tus by splitting the OptS and SubS into two: grouped and 
ungrouped galaxies. We then compare the two samples and 
find that the SubS appears to have galaxies which are al- 
most equally split between the field (52%) and groups (48%), 
whereas the OptS shows more of a difference with 62% in 
the field and 38% in groups. ETGs within the groups can 
be further split according to the multiplicity of the group: 
ETGs in small groups with less than five galaxies (32% and 
24% for the SubS and OptS, respectively), and larger groups 
with five or more galaxies (15% and 14% respectively). 

Fig. [11] shows the distribution of surface densities with 
stellar mass for our two samples. The galaxies have been 
coloured by the multiplicities of their groups, making it sim- 
ple to identify lone ETGs, ETGs in small groups and those in 
large groups. Both plots show that higher densities are linked 
with larger groups, but are not exclusively dominated by 
these larger groups. Both of these plots also show there are 
some ETGs with high surface densities which do not seem to 
share a group. This may be an effect of the groups data selec- 
tion function, which changes with redshift. This figure high- 
lights the similar SubS and OptS mass ranges, but different 
galaxy surface density ranges. Testing for Spearman corre- 
lations does not reveal any trends between surface density, 
group multiplicity and stellar mass, although these same 
tests suggest that the OptS contains more ungrouped ETGs 
in high density areas (rs=0.005) than the SubS (rs=0.24). 

There also appear to be some missing galaxies in Fig. 
Illf a). The OptS shows galaxies of quite low stellar mass 
(M* <1O'^°M0) going to higher densities, but these do not 
appear for the SubS. In general, H- ATLAS is known to pref- 
erentially detect higher mass galaxies, which could partially 
explain this. In addition, the larger H-ATLAS PSF when 
compared to that of the optical may result in low coun- 
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Figure 11. Variation of stellar mass with surface density, (a) 
shows results for SubS and (b) for OptS. Both plots are coloured 
by group multiplicity: open blue circles are galaxies which are un- 
grouped and assigned a group size of 1, red filled circles represent 
ETGs in small groups (< 5 galaxies) and large, green filled cir- 
cles those in large groups (> 5 galaxies). Note that the horizontal 
scales differ in these two panels. 



terpart reliability for ETGs in the densest regions, result- 
ing in their removal from our samples. The sparsity of low 
mass, dusty ETGs at high surface densities may be caused 
by sample completeness effects. However, the sparsity of 
dusty ETGs in general at high densities suggests that dusty 
ETGs do not occupy the densest regions o f groups and clus- 
ters, u nlike their non-dusty counterparts. iKauffmann et all 
l|2004 ) showed that generally, massive galaxies in low den- 
sity environments in their SDSS sample con tain the most 
dust, based on evidence of optical attenuation. iKavirai et al\ 
studied optical images of ETGs with visible dust 
lanes and patches. They found that dusty ETGs occupy less 
dense environments than those with no sign of dust obscu- 
ration. Our findings are qualitatively consistent with these 
results. 



5 DUST PROPERTIES OF DETECTED ETGS 

5.1 Fitting Modified Planck Functions 

The SubS ETGs have both Herschel FAGS 100 and 160^m 
measurements, as well as 250, 350 and 500pim fiuxes from 
the SPIRE instrument. We select a subsample (3a sample) 
of 188 ETGs which have at least 3cr flux levels (22.6 mjy) 
in the 350/im SPIRE waveband, as well as the requisite 5a 
emission in at least one SPIRE sub-mm waveband (see i]3.4p . 
PACS measurements are missing for 27 of these galaxies, 
and for these we choose to fit only SPIRE flux densities. As 
specified for all H- ATLAS Phase 1 data, fiux errors of 7% of 
the catalogue fiux values are added in quadrature to account 
for the uncertainty in the SPIRE photometric calibration, 
and 10% to the PACS fluxes. 

In order to estimate dust temperature and mass we then 
fit isothermal, modified Planck functions to the data for each 
galaxy, of the form: 



where Bx is the Planck function, Q is the amplitude param- 
eter for the model fit, a nd P the dust emissivity index which 



we fixed as 2.0 (e . g. Galametz et al.l I2OIII : iDavies et al.l 



l2012l : ICortese et all 120121 ). This grevbodv model is fit to 
the FIR/sub-mm fiux densities, which are corrected for red- 
shift. We are thereby able to measure the typical rest-frame 
temperatures of emission by the dust distributions in these 
galaxies. 

Dust masses for this subsample are calculated using 



F250DIK 



«250B(T)250 (1 + z) 



(10) 



(IWhittet|[l993 . Eq. 6.12), where F250 is the observed fiux 
value in the 250/im waveband in Janskies (Jy), Dl the lu- 
minosity distance to th e source and K is t he K-correction 
as defined in Eq. 2 of iDunne et al.l (|201ll ). The assumed 
mass absorption co efficient (K250) is 0.89m^kg~^ at 250/im 
dDunne et al.ll201ll ). The value of B(T)250 in Eq.[lO]is cal- 
culated using the best fit, cold dust, rest-frame temperature 
from the SED fitting for each galaxy. We choose to measure 
the dust mass using this waveband because of the relatively 
smaller errors when compared with the other SPIRE wave- 
bands. 

The results of our fitting give the mean cold dust tem- 
perature for the ETGs as 22.ll|^K with 5*'' to 95*'' per- 
centiles of 14-28K, with no obvious variation of fit tem- 
perature with stellar mass. The mean cold dust mass is 
calculated as (1.8±0.5) x lO'^Mo with a percentile range 
of (0.19-5.41)x10''Mq. The resultant average specific dust 
mass (logio(Md/Mt)) is calculated as -3.37. Errors on tem- 
perature are given as formal uncertainties and propagated 
through Eq. [TO] Due to the non-linearity of this propaga- 
tion, we only show average dust mass errors for those galax- 
ies with fuU PACS and SPIRE data. Thus the correlation 
between dust temperature and mass is accounted for within 
dust mass errors. 

We present our specific dust masses in Fig. 1121 plot- 
ted against stellar mass. Overall, the galaxies show a clear 
trend for lower mass galaxies to contain higher normalised 
dust masses, with a straight line fit to the distribution yield- 
ing a Pearson coefficient (r^) of -0.55. The dependancy on 
morphology within ETGs is addressed in Here we com- 
pare our average normalised dust mass with those from other 
works with Herschel. These results and the standard errors 
on the mean are also shown in Fig. I12f right) for ellipticals 
and lenticulars separately. 

Our ETC sample does not contain many low mass 
ETGs, which indicates that there may be missing galaxies 
at the low stellar mass and specific dust mass end of Fig. 1121 
To test for this, we show a cross-hatched region on the plot, 
representing the 250/im 5a limit. This region is based on the 
same dust mass calculations as our ETGs, with a represen- 
tative upper limit dust temperature of 30 K and minimum 
redshift of 0.013. This excluded region becomes larger with 
decreasing temperature and increasing redshift, and repre- 
sents the region of the plot which H- ATLAS is insensitive to. 
This indicates that we are unlikely to be missing ETGs on 
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Figure 12. Specific dust mass of ETGs in our sample SubS as a function of stellar mass. The points are coloured by morphology: dark 
blue filled circles for ellipticals (E) and light blue open circles for lenticulars (SOa). Average errors are shown as bars on the top right. 
Same colour dashed lines represent best fit straight lines through the respective points. The black solid line is the overall fit to all points. 
The 250/xm flux limit for a typical galaxy with a upper limit dust temperature of 30K is indicated by the grey cross-hatched region. 
In the panel to the right, filled diamonds and error bars represent mean specific dust masses for a range of studies with Herschel data, 
described in the text. 



this diagram due to the flux hmit, since we would expect to 
detect an increasing proportion of galaxies above this limit. 

Based on the results found in i|4j we want to further 
explore the properties of the ETGs for high and low lumi- 
nosity subsets. Luminosity and stellar mass are known to be 
linked, and we plot these properties to find that values of 
Mr=-21.5 mag / M*=1O^°-^M0 can be used as a high and 
low luminosity /stellar mass divider. From Fig. [12] we see 
that while lenticulars dominate the median range of stel- 
lar masses, ellipticals appear to be separated into low/high 
stellar masses. Those ETGs with the highest stellar masses 
show the lowest specific dust masses. We have shown that 
the brightest ETGs contain Sersic indices most similar to 
undetected ETGs (see Q and this therefore indicates that 
ETGs with low specific dust masses have high Sersic indices 
and are not dissimilar to ETGs in the OptS in this respect. 

We next compare results with ot her H-ATLAS da t a. Re - 
cent multi-wavelength SED fitting bv lRowlands et al.l (I2OI2I ') 
on a sample of 42 detected ETGs in the H-ATLAS SDP 
gave average cold dust masses as S.SxlO^M© for an al- 



lowed dis tribution of temperatur es between 15-25K. Improv- 
mg upon iRowlands eTsI} (|2012l '). we have chosen to fit the 
ETGs with 3cr detections at 350/^m, and we have chosen 
a more straightforward single-component fitting approach. 
This low-redshift SubS alone is five times the number of 
galaxies in the SDP and therefore a more statistically sig- 
nificant result. Our resul ts indicate lower mean specific dust 
masses for ETGs than in lRowlands et sH (|2012l ). which may 
be partly due to the nearer redshift limit of our SubS sam- 
ple presented here. For the sake of comparison, we have also 
shown the average normalised dust mass for their larger sam- 
ple of LTGs. This is clearly a lot higher than the ETGs, as is 
to be expected for spiral galaxies, which are known to have 
high dust masses. 

We can also compare our results w i th tho se from the 
KINGFISH survey, where ISkibba et all (|201ll ) also calcu- 
late dust temperatures and masses for a sample of 10 nearby 
ETGs using single-temperature modified blackbodies. Their 
specific dust masses are lower than ours, but fall well within 
the range of masses exhibited by our galaxies. Therefore we 
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Figure 13. Distribution of specific dust masses with (NUV-r) 
colour for tlie SubS. Points are coloured by the log of their stellar 
mass. The solid black line is the best fit straight line through the 
points, with the Pearson coefficient for the fit labelled in the top 
right-hand corner. Typical error bars for the points are shown in 
the top left corner. 



can attribute the difference in mean masses to sam ple size. 
This is less of an issue with the lSmith et al.l (|2012l ) study in 
the HRS, where they select a sample of 62 ETGs, mainly in 
the Virgo Cluster, of which 45% are detected in the sub-mm. 
They also calculate dust masses for the 28 detected ETGs 
using single-temperature greybody fitting, and their average 
specific dust mass for detected ETGs is much lower than any 
of the aforementioned works, including ours. However, they 
do not seem to sample as large a range of environments as 
these other works, since their sample is restricted mainly to 
galaxies in the high densities of the Virgo Cluster. Further 
HRS work by ICortese et all l|2012l ) shows that in clusters 
such as Virgo, the dust fraction of galaxies of a given mass 
is significantly lower than that of galaxies in isolation or in 
small groups. This bias for dusty ETGs to reside in lower 
density environments is also hinted at in Table [5] and Fig. 
[11] for our SubS. The KINGFISH and HRS studies look 
at galaxies in the very nearby universe, and therefore they 
are able to detect the lowes t dust masses in galaxies. Addi- 
tionally, ISmith et al] (|2012l ) emphasise that they only study 
the most massive ellipticals in the nearby universe. Given 
this information, the differences in average normalised dust 
masses between this and these other two Herschel works 
appear quite small. 

Next we consider how the computed dust masses vary 
with UV-optical colours. This is presented in Fig. 1131 where 
the galaxies are also coloured by stellar masses. Here we 
identify a very clear trend. As our galaxies get bluer, they 
also increase their specific dust masses and decrease in stellar 
mass. This highlights a key result: the most massive ETGs 
are redder and therefore more quiescent, and they also con- 
tain proportionally less dust. This strengthens the link be- 
tween dust and star formation and, assuming the more qui- 
escent ETGs do not host many young blue stars, predicts 



that low mass red stars do not contribute as highly towards 
the dust presence in ETGs as brig ht young stars . This sort of 
evolution has been shown before: iBourne et all ||2012D show 
a similar strong trend for their sample of red galaxies. To 
further understand it, we would have to work out whether 
we are looking at these ETGs at different ages, and therefore 
different stages of their evolution, or whether we are looking 
at two completely different populations and types of ETGs. 

5.2 Contamination Issues 

iRigbv et al.l (|201ll ) found that a significant number of 
sources may have 350 and 50Q/.im flux densities that are over- 
estimated by a factor of ~2 in the H- ATLAS SDP fields. This 
is due mainly to source confusion where the signal-to-noise 
levels are low. This would result in SED fitting underesti- 
mating dust temperatures and overestimating dust masses. 
There is nothing that can be done to remedy this on an 
object-by-object basis, but we should be aware of such an 
effect having carried through to the Phase 1 data. 

We consider whether any of our ETGs are radio- 
emitting AGN by referring to the Very Large Array FIRST 
survey, which cover 10,000 degrees over the North and South 
Galactic Ca ps, currently with 30% optical counterparts in 
the SDSS (|Becker et al.1 119941 ) and covering the GAMA 
equatorial regions down to the SDSS DR6 brightness limit. 
We searched for potential counterparts within a 20" search 
radius and found two of our SubS ETGs, within search radii 
of 18.9" and 9.4", and 1.4 GHz integrated fiux densities of 
1.17 and 3.59 mjy respectively. Given the likelihood of these 
being false counterparts and the low associated flux densi- 
ties, we deem radio emitting AGN to not be an issue for our 
ETGs. 

If the ETGs do host radio-emitting AGN which we have 
been unable to detect through FIRST counterpart matching, 
it is necessary to consider their implications in our results. 
Synchrotron emission from the AGN would be modelled as a 
radio power-law, which may extend into the longest sub-mm 
wavebands. Our greybody fit parameters will then represent 
not only the thermal FIR/sub-mm emission, but also the 
synchrotron emission component. This could lead to false 
cooler temperatures resulting from our fits, as the apparent 
greybody peak would be pushed to longer wavelengths. 

We double-check whether radio-loud AGN are present 
in our ETGs by looking at t he SPIRE 250/350 and 350/500 
colours. Bosel li et all l|2010l ) show that SPIRE colours are 
useful for discriminating thermal from synchrotron emission 
in radio galaxies, with the colours becoming very small for 
those galaxies with synchrotron emission. The low redshift 
sample we are fitting in this section have high (f350/f500>l) 
SPIRE ratios, providing further evidence that the 500/im 
waveband is unaffected by synchrotron emission in our sam- 
ple. 

Free-free emission (thermal bremsstrahlung) from 
ionised HII regions is also a potential contamination fac- 
tor in the 500/im waveband. However, because ETGs are 
mostly quiescent, such emission from these regions is likely 
to be small and is difficult to quantify; 

As discussed in |]3.4I w e used the iNegrello et all (|2010l ) 
and IConz alcz-Nuc vo et al.l (|201^ selection criteria to re- 
move nine galaxies which are quite likely to be lensed. These 
galaxies did have high enough flux emission to run grey- 
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body fits across them. The results of these fits were a range 
of dust masses between 10''"*Mq (-3.1<logio(Md/M.) <- 
1.9) and temperatures between 20-35K, with only one lens 
galaxy candidate given a temperature as low as IIK. Al- 
though these galaxies are not included within our dust mass 
plots, upon inspection we note that if they were included 
they would not change the straight line fits in Figs. 1121 and 
1 131 or affect the results significantly. 

Finally, as a check, we re-ran the greybody fitting on 
the ETGs for just the four shortest waveband fluxes, thereby 
excluding the 500/im points. The resultant rest-frame tem- 
peratures did not change greatly for many of the galaxies, 
except in the cases (~4%) where flux from the 250/im band 
was much greater than the 350/im flux. The lack of a strong 
systematic bias occurring when running this test gives fur- 
ther evidence that the dust characteristics of our SubS are 
not significantly affected by synchrotron emission or other 
contamination issues affecting the 500/^m band fiuxes. 

5.3 Elliptical vs Lenticular Dust Characteristics 

It is necessary to consider the implications that SOa galaxies 
have created by being included within our samples, as they 
may have had an effect on dust results. 

It is possible that lenticular galaxies are in fact Spi- 
rals which have had their s tar-formation cut off (e.g. 
lAragon-Salamanca et all |2006l). Supposing ellipticals are 
formed in a different way, possibly via mergers, the two dif- 
ferent formation methods will result in different proportions 
of dust left within the galaxy systems. If this is the case, we 
are left with one of three scenarios. Either the dust masses 
and temperatures we are fitting to our sample are repre- 
sentive of dusty ETGs overall, or they are representative of 
either dusty SOa galaxies or E galaxies separately. 

Fig. [12] shows the trend of normalised dust mass with 
stellar mass for all ETGs, coloured by morphology. We plot 
these separately because lenticular and elliptical galaxies are 
not necessarily at the same stage of their lives, and may have 
different evolutionary patterns. Additionally, the presence 
of a disk in lenticular galaxies points towards them possibly 
having higher dust masses than ellipticals, thus skewing our 
perception of the dust contents of elliptical galaxies. Such a 
result is not indicated in Fig. 1121 where the lenticulars are 
shown to have the same mean specific dust mass (-3.35) as 
the ellipticals (-3.36). However, a KS-test of the stellar mass 
distributions show them to be marginally different at the 1% 
level (KS-prob=0.009). The straight-line fit to the ellipticals 
suggests a slightly stronger trend for specific dust mass to 
change with stellar mass (rp=-0.71) than the trend for lentic- 
ulars (rp=-0.41). Smaller number densities for the ellipticals 
weakens this trend somewhat. This result nevertheless indi- 
cates that the two different types of ETGs have similar dust 
properties and therefore that there is some similarity in the 
way dust in these galaxies evolves and consequently, in the 
evolution of the galaxy itself. 

We explore the separate morphological properties fur- 
ther in Figs. [T?] and [TJ] by plotting different parameters 
against the absolute r band magnitudes and normalised dust 
masses respectively. In Figs. ll4f a') and (b) we show how both 
normalised dust mass and total dust mass vary with abso- 
lute magnitude. For both galaxy types we see a clear trend 
emerging here, which is similar to the trend identified earlier 



in ^5.11 This is expected since magnitude and stellar mass 
are related (see discussion in ^STTJ. These plots show that 
the brightest ETGs contain proportionally the least dust 
and vice versa (Fig. I14f d')). Once again, by breaking ETGs 
down to elliptical and lenticular, we explore their compar- 
ative properties and find the brightest ETGs are elliptical, 
with low specific dust mass (Fig. [Mjd)) and high dust mass 
(Fig. I14f e')). On the faint end, we see a well mixed distribu- 
tion of both types of ETGs. This, and a KS-test, indicates 
no preference for lenticulars to have different specific dust 
masses than ellipticals. 

The wide distribution of dust temperatures arising from 
our modified Planck function fits is also highlighted in Fig. 
I14r c) and (f). Our ETGs extend to very low dust tempera- 
tures, with errors of typically ±2 K. This could be attributed 
to the fact that this is a very large sample of ETGs observed 
by very sensitive detectors and the galaxies may contain dust 
temperatures which could not physically have been mea- 
sured up until now. However, it should be noted that the 
galaxy with the lowest temperature (~9 K) has large PACS 
measurement errors. 

All the other ETGs in our sample have dust tempera- 
tures higher than 10 K. The distributions of ellipticals ver- 
sus lenticulars look quite similar, although the lowest tem- 
peratures (<15 K) belong predominantly to ellipticals. KS- 
testing shows the distributions are not significantly different. 
Based on this result, it appears that lenticulars not only 
contain similar (specific) dust masses to ellipticals, but also 
similar dust temperatures. The lack of observable trend of 
temperature with luminosity also implies that there is no di- 
rect correlation between distributions of dust temperature 
and stellar mass for the two types of ETGs. 

The distributions of environmental, Sersic and colour 
properties of the galaxies as a function of normalised dust 
mass are examined in Fig. 1151 Fig. I15f a) indicates that there 
is no significant difference in the distribution of surface den- 
sities for the two morphologies, which is confirmed by a KS- 
test. Its counterpart scatter plot (Fig. llSf d)) shows no ob- 
servable trend for these properties, although the few ellipti- 
cals with an intermediate range of specific dust mass seem 
to be clustered towards the lower densities. 

Figs. I15f b) and (e) show the Sersic (ur) distributions 
of these morphologies and highlight a useful characteristic 
of ellipticals. Although there are only small numbers of el- 
lipticals, approximately a third of them are at lower Sersic 
(n~l). However, Fig. llSf el also indicates that on average, 
these low Sersic ETGs are those with the highest specific 
dust masses (rp=0.1). This again supports the idea that 
du st is causing the low ering of Sersic index, as predicted 
bv lPastrav et all l|2013l ). 

Finally we examine UV-optical colour in Figs. llSf c'l and 
(f). Although this trend was seen previously in Fig. 1131 in 
this plot we see the significance of different morphologies. 
Once again we find different groupings of ellipticals, unlike 
the solidly filled range presented by the lenticulars. A KS- 
test shows NUV-r distributions to be significantly different 
at well below the 1% level (KS-prob=0.0009). Elfipticals are 
grouped either as very blue and dusty, or red with a range 
of specific dust masses. Such a distribution of specific dust 
mass at the red end indicates that dust may be contributing 
to reddening in this plot, otherwise we would expect galaxies 
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Figure 14. Scatter plots and normalised histograms for SubS ellipticals and lenticulars. Ellipticals are represented by dark blue filled 
circles and filled histograms, and lenticulars by light blue filled circles and histograms. From left to right: (a) distributions of dust mass 
normalised by stellar mass, (b) histograms of logged dust mass, (c) histograms of rest-frame dust temperatures calculated by modified 
Planck function fitting, (d) absolute r band magnitude (Mr) against normalised dust mass, (e) M,- against logged dust mass, and (f) M,- 
plotted against best fit dust temperature. Average errors are shown in the top-right corners. KS-probabilities described in the main text 
are also included for each histogram. 



with such high levels of dust to occupy the bluer end of this 
distribution. 



6 SUMMARY AND CONCLUSIONS 

Tliis work is the first in a two-part series exploring intrinsic 
and dust properties of early-type galaxies. By using visual 
classi fication s applied to a low-z sample of GAMA galaxies 
from iKelvin et al.l (|in prepi) , we created two ETG samples 
from the GAMA data. These are the SubS - galaxies clas- 
sified as ETGs with H-ATLAS SPIRE 5a detections - and 
the Opts - galaxies classified as ETGs without H-ATLAS 
SPIRE detections. 

We explored a series of optical and UV properties for 
these two samples. These included concentration and Sersic 
index, stellar masses, luminosities, optical colours, and ef- 
fective radii. GALEX NUV detections were used to explore 
UV-optical colours of our galaxies and the environments of 
the ETGs were examined by way of 5"' nearest neighbour 



densities and group membership data. The results for these 
tests are summarised here: 

i) We find a 29% H-ATLAS detection rate for ETGs in the 
equatorial field. The detection rate for pure sub-mm de- 
tected ellipticals is 10% (see i]3.3|) . 

ii) H-ATLAS detected {SubS) ETGs are shown to have lower 
concentrations and Sersic indices than undetected (OptS) 
ETGs. It has been shown that th e presence of dust can lower 
Sersic index ( jPastrav et al.ll| 2013i ), an effect that we consider 
to be a possible explanation for our findings. This effect 
will be further explored in Paper II of this series, where we 
will investigate whether this is true for our samples in all 
wavebands and at higher redshifts. 

iii) H-ATLAS detected ETGs are shown to be typically op- 
tically brighter, bluer and larger than undetected ETGs. 
However, differences in the stellar mass distributions of the 
samples are found to be not significant. 

iv) Our H-ATLAS detected sample shows significantly brighter 
NUV luminosities and bluer colours in the UV-optical 
than the undetected sample. This colour trend is in agree- 
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Figure 15. Scatter plots and normalised histograms for SubS ellipticals and lenticulars. Ellipticals are represented by dark blue filled 
circles and filled histograms, and lenticulars by light blue filled circles and histograms. From left to right: (a) normalised logged distri- 
butions of surface density E;;, (b) logged histogram of r band Sersic index, (c) distributions of UV-optical colour (NUV-r), (d) scatter 
plot showing specific dust mass against surface density, (e) specific dust mass as a function of Sersic index, and (f) as a function of 
NUV-r colour. Average errors are shown in the top-right corners. KS-probabilities described in the main text are also included for each 
histogram. 



me nt with previous wor ks such as [Rowlands erall (|2012l ) 
and lDariush et all (|201ll '). However, the H- ATLAS detected 
ETGs are redder than LTGs from the same parent sample, 
and therefore dominate an intermediate NUV-r colour space, 
v) The H- ATLAS detected ETGs are shown to inhabit sparser 
environments, particularly for the fainter ETG subset. By 
comparing our ETGs in both small and large galaxy groups 
with ungrouped ETGs, we find that for the SubS, higher 
mass stellar systems are typically associated with larger 
group sizes and denser environments. 

To investigate dust properties, we fitted modified single- 
temperature Planck functions to the 250/im 5a detected 
SubS galaxies containing >3a emission in the SPIRE 350/im 
waveband. We then presented characteristics of the SubS 
based on resultant best-fit temperatures and associated dust 
masses, and explored trends with optical and UV properties. 
These results can be summarised as follows: 

i) ETGs in our H- ATLAS detected sample are shown to con- 
tain dust masses ranging from 8. 1x10^-3. 5xlO*M0, with a 



range of rest-frame temperatures from 9-30K. These dust 
masses are consiste nt with previous Herschel work such as 
ISkibba et al.l (1201 ih. but lower than resu lts from previous H- 
ATLAS work in [Rowlands et"al1 (120121 '). These resuhs may 
differ due to our larger sample size and our lower redshift 
limit. 

ii) We discover a strong trend for specific dust mass to de- 
crease with redder (NUV-r) colour. This implies that the 
dustiest ETGs have the bluest colours, linking recent star 
formation with a heavier dust presence. 

iii) We show that the faintest H- ATLAS detected ETGs with 
the lowest stellar masses have higher specific dust masses 
and lower Sersic indices compared to the brightest H- ATLAS 
detected ETGs. This result is consistent with downsizing, as 
the most massive, brightest galaxies are more similar to the 
Opts ETGs in terms of their specific dust masses and Sersic 
profiles. 

iv) By splitting our dusty ETG sample into separate morpholo- 
gies of elliptical and lenticular galaxies, we gauge how they 
contribute to the measured dust properties. We find no dif- 
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ference in the dust masses (whole or specific) or tempera- 
tures of both ETG types. ISowever, we find significant differ- 
ences in their NUV-r colour and stellar mass distributions, 
v) We show that ellipticals may be grouped into two sets: the 
faint, blue, low mass, relatively dusty ellipticals and bright, 
red, massive, ellipticals with lower specific dust masses. This 
leads us to the conclusion that we may be studying two 
different populations of ellipticals, or at least two different 
age ranges. 

We intend to further explore these separate populations 
of ETGs by pushing our samples out to higher redshifts in 
Paper II of this series. As we can no longer rely on visual 
classifications at these redshifts, we will explore proxies for 
morphology. We will then repeat our study of ETG optical 
and dust characteristics, and compare them with these low 
redshift samples. 
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